Spin and temperature dependent study of exchange and correlation in thick 

two-dimensional electron layers. 
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The exchange and correlation E xc of strongly correlated electrons in 2D layers of finite width are 
studied as a function of the density parameter r 3 , spin-polarization £ and the temperature T. We 
explicitly treat strong-correlation effects via pair-distribution functions, and introduce an equivalent 
constant-density approximation (CDA) applicable to all the inhomogeneous densities encountered 
here. The width w defined via the CDA provides a length scale defining the z-extension of the quasi- 
2D layer resident in the x-y plane. The correlation energy E c of the quasi- 2D system is presented as 
an interpolation between a ID gas of electron-rods (for w/r s > 1) coupled via a log(r) interaction, 
and a 3D Coulomb fluid closely approximated from the known three-dimensional correlation energy 
when w/r 3 is small. Results for the E xc (r s ,C,,T), the transition to a spin-polarized phase, the 
effective mass to*, the Lande g-factor etc., are reported here. 

PACS numbers: PACS Numbers: 05.30.Fk, 71.10.+a, 71.45.Gm 



INTRODUCTION. 



The 2D electron systems (2DES) present in GaAs or 
Si/SiC>2 structures access a wide range of electron densi- 
ties, providing a wealth of experimental observations 0. 
The 2D electrons reside in the x-y plane and also have a 
transverse (^-dependent) density n[z) which is confined 
to the lowest subband of the hetero-structure 2] . The 
2D character arises since the higher subbands are suffi- 
ciently above the Fermi energy Ep. Then the physics 
depends only on the "coupling parameter" T = (poten- 
tial energy)/ (kinetic energy), the z-distribution n(z), the 
spin-polarization £, and the temperature T which has 
to be significantly smaller than the Fermi energy Ep to 
preserve the 2-D character. The T for the 2DES at the 
density n is equal to the mean-disk radius r s = (7m) _1 / 2 
per electron, expressed in effective atomic units which 
depend on the bandstructure mass nib and the "back- 
ground" dielectric constant t\>. The z-motion in the 
lowest subband may have widths of ~ 600 A, in GaAs 
when r s is ~ 6, in heterojunction insulated-gate field- 
effect transistors (HIGFET) which have been an object 
of recent studies Q]- Similarly, other nanostructures (e.g., 
quantum dots) contain electrons confined to a micro- 
region in the x-y plane, and have a sizable z extension^]. 
Hence layer-thickness effects are important in many ar- 
eas of nanostructure physics. Appropriate correlation 
functionals|5j for such systems are still unavailable, even 
though the exchange functional for Fang-Howard distri- 
butions is known |(j . 

Layer thickness effects are a long standing probe of ex- 
change and correlation theories in 3D electron slabs[7j]. 
The relevance of the finite size of the 2D layers had also 
been considered within the quantum Hall effect HE|, and 



more recently in the context of the q-fa ctor and the ef- 
fective mass m* of 2-D layer ElM El El E3 In 
the early days of the application of diagrammatic per- 
turbation theory (PT) to condensed-matter problems, it 
was normal to attempt to calculate various many-body 
properties like the effective mass m* , the effective g- 
factor g*, and corrections to the total energy using per- 
turbation methods. The need to go beyond the random- 
phase approximation (RPA) was rapidly appreciated and 
lead to the work of Hubbard, Rice, Vosko, Geldart and 
others El- The common experience with the generalized 
RPA method, when applied to 3D electrons and to ideally 
thin 2D layer is that it predicts a spin-polarized phase at 
unrealistic high densities (low r s ), while Quantum Monte 
Carlo (QMC) simulations suggest a density near r s ~ 25- 
26 in 2D. RPA methods lead to negative pair-distribution 
functions (PDFs), incorrect local-field corrections in the 
response functions, and disagreement with the compress- 
ibility sum rule and other formal conditions. Recent 
attempts to calculate the effective mass m* for ideally 
thin layers El Ei 

show strong disagreement with the m* 
obtained from QMC simulations 17] . In fact, the main 
thrust of the programs of Singwi, Tosi et al. (STLS)El, 
and Ichimaru et al. E3 was to overcome the shortcomings 
of the RPA-like approach via non-perturbative methods. 
Many-body calculations where different parts of the cal- 
culation (e.g, vertex corrections, local-field corrections, 
etc.) are obtained from different sources, (e.g, vertex 
corrections from a model, local-field corrections or cor- 
relation energies from QMC, and some other quantities 
fitted to sum rules etc.) and combined together, have 
also appeared. Unless the same "mixture" is used to cal- 
culate a multitude of properties and shown to lead to 
consistent results, such theories have to be treated with 
caution. 
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We have introduced an approximate method for 
strongly correlated quantum systems where the objec- 
tive is to work with the PDF of the quantum fluid, 
generated from an equivalent classical Coulomb fluid 
whose temperature T q is chosen to reproduce the cor- 
relation energy of the original quantum fluid at T = 0. 
The classical PDFs are calculated via the hyper-netted- 
chain (HNC) equation, and the method is called the 
CHNC. As the method has been described in previous 
publications [jfl Ell I22 I Ei l . and successfully applied to a 
variety of problems |24l l25l Erj only a brief account 
is given here, mainly to help the reader. The PDFs are 
obtained from an integral equation which can be recast 
into a classical Kohn-Sham form where the correlation ef- 
fects are captured as a sum of HNC diagrams and bridge 
diagrams. 

-/9(P«(r)+V eo „(r)+JVy(r)+B«(r)) 



gij(r) = exp- 



(1) 



The temperature of the classical fluid T q = 1/(3 is chosen 
such that at T = the calculated classical g(r) recovers 
the known correlation energy of the fully spin-polarized 
2D electron fluid at the given density. This fitting has 
been done in ref. , and T q is known as a function of r s . 
At finite-T, the classical fluid temperature T c f is taken 
to be 



T cf = (T ? 2 + T 2 ) 1/2 (2) 

This T c f is used for all spin polarizations. In Eq. ^the 
indices i,j label the spin states. The Pij(r) is chosen to 
ensure that gtj(r) reduces to the explicitly known non- 
interacting PDF, gfj(r), when the Coulomb interaction 
V cou (r) is switched off. Thus Pij(r) takes care of Pauli- 
exclusion effects and ensures that the "Fermi-hole" is ex- 
actly recovered [isj. Also, Nij(r) is a sum of terms that 
appear in the hyper- netted-chain equation, while Bij(r) 
contains 3-body and higher-order diagrams not contained 
in the nodal term Nij(r). The latter depends implicitly 
on gij(r), and is evaluated via the Ornstein-Zernike equa- 
tion. The bridge term is very difficult to evaluate, but the 
hard-sphere fluid provides a good approximation. That 
is, in 2D, we specify Bij(r) by specifying the hard-disk 
radius rjy, or equivalently the packing fraction 77, and 
use the Percus-Yevik approach|29j. As the parallel-spin 
three-body clusters are suppressed by the Pauli exclusion, 
we use only a single bridge function, viz., Bi2(r). This 
makes the bridge interaction effectively independent of C- 
Bulutay and Tanatar, and also Khanh and Totsuji30], 
have examined variants of CHNC without a bridge func- 
tion (this is some what equivalent to neglecting back-flow 
terms in QMC simulations of 2D systems). The hard- 
sphere radius rn, and the packing fraction 77 are given 
by|12: 



r H = 2r s i 1 1 l 2 

r, = 0.1175ry 3 (t c/ +t 2 /2) 
t c f = T c f/Ep, t = T/Ep 



(3) 
(4) 
(5) 
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FIG. 1: The bridge function B{r) for r s — 3 and r s 
effect of changing the temperature is also shown. 



5. The 



A plot of the bridge function for a few typical cases is 
given in Fig. ^ 

The CHNC method was applied to the 3D and 2D elec- 
tron fluidsHE IH |H i3, to dense hydrogen fluidQ, 
and also to the two-valley system in Si/Si02 2DES;25, 
Et1 |. In each case we showed that the PDFs, energies and 
other properties obtained from CHNC were in excellent 
agreement with comparable QMC results. The advantage 
of the CHNC method is that it affords a simple, semi- 
analytic theory for strongly correlated systems where 
QMC becomes prohibitive or technically impossible to 
carry out. The classical- fluid model allows for physically 
motivated treatments of complex issues like three-body 
clustering etc., which are difficult within quantum meth- 
ods. The disadvantage, typical of such many-body ap- 
proaches, is that at present it remains an "extrapola- 
tion" taking off from the results of a model fluid. The 
fully spin-polarized infinitely-thin uniform 2DES is the 
model fluid[i3|, as the E c of this one-component system 
is accurately known. 

In this study we use a method of replacing the inhomo- 
geneous electron distribution n(f , z) = n(r)n(z) by a ho- 
mogeneous density system (l5H32T | i.e., a constant-density 
approximation (CDA) where the transverse density n(z) 
is a constant within a slab of width w, and zero out- 
side. The CDA avoids difficulties associated with gradi- 
ent expansions noted in ref.[5j- Also, the CDA presents a 
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unified approach to quasi-2D distributions like the Fang- 
Howard model0, the quantum-well model etc. E xc for 
such distributions is calculated as a function of the spin 
polarization £, 2-D density parameter r s , the CDA width 
w and the temperature T. This enables us to determine 
physical quantities related to the Landau Fermi liquid 
parameters. Thus the spin susceptibility enhancement, 
the effective mass m* , and the Lande g factor for the 
quasi 2D electrons are presented. 



THE QUASI-2D INTERACTION 

The transverse distribution n(z) of the 2D electrons is 
given by the square of the lowest subband wavefunction 
4>{z) of the heterostructure, calculated within the enve- 
lope approximation. The nature of the materials used 
(e.g, Si/SiC>2 or GaAs) and the doping profiles determine 
the confining potential and the electron density n(z) in 
the z-direction. Typically, n[z) may be modeled by one 
of the following forms. 



n(z) — S(z), ideal thin layer 

= (2/w) sin 2 (7rz/w), infinite well 
= (b 3 /2)z 2 exp(— bz), Fang-Howard 



(6) 
(7) 
(8) 

= 1/w, \z\ < w constant-density model (9) 

The second and third are frequently used approximate 
(but adequate) models, while we present the fourth 
model, the CDM. This is an excellent constant-density 
approximation (to be called the constant-density approx- 
imation, CDA) to generate the effective 2-D potential 
arising from most n(z) distributions, on suitably choosing 
w. The Fang-Howard (FH) form|2L l31|. contains the pa- 
rameter 6, and is normalized in the range to oo. The FH 
parameter b is such that b 3 — (48 * ir/a )(rid + lln s /32), 
where n s is the 2D electron density n, and rid is the de- 
pletion charge density 0). The material parameters are 
contained in the effective Bohr radius clq = ebh 2 /rribe 2 
defined in terms of the usual constants, e& and mb being 
the dielectric constant and the band mass. For the de- 
vices used in ref. 0, [n) , the depletion density has been 
reported to be negligible 12). Then b 3 = 33/(2r 2 ), in 
atomic units. 



The constant-density model. 

We denote the Coulomb potential in an infinitely thin 
layer by V(r) = 1/r, while W(r) is used for the effective 
2-D potential of a thick layer. The effective 2D-Coulomb 
potential Wir) between two electrons having coordinates 
{fx, z\) and (fi, z%), with r = rj — r5 is given by, 



Here z m is 00 for FH, while z m — w for the others. The 
potential W(r) — (l/r)F(r) and the form factor F(r) 
reflects the effect of the z-extension of the density. There 
is no analytic form for F(r) in the Fang-Howard case, 
while the g-space form, F(q) is knownQ- If the dielectric 
constants of the barrier and well material were assumed 
equal, then the Fang-Howard form factor is: 

Here we derive a potential W(r, w) for the constant- 
density model(CDM) which is, to an excellent approx- 
imation electrostatically equivalent to the the 2D poten- 
tial for any reasonable n(z). These FH-type distribu- 
tions are themselves convenient fits to the self-consistent 
Schrodinger solutions and have uncertainties of a few per- 
cent. The CDA holds well within such limits. The poten- 
tials are explicitly shown in Fig.[21for the FH form where 
we have taken an extreme example with b = 0.1. The 
method of replacing an inhomogcncous distribution by a 
uniform distribution is suggested by the observation that 
the non- interacting total pair-correlation function h°(r) 
has the form ~ n(r) 2 , where n(r) is the density -profile 
around the Fermi hole. In our case we wish to replace 
the inhomogeneous n{z) by a constant-density distribu- 
tion n c d which has the same electrostatic potential in the 
2-D plane as n(z). 



n c d = 1/w = / n(z) 2 dz 



(12) 



Since the subband distribution is normalized to unity, 
the width w of the constant-density slab is simply 1 / n r A. 
Starting from different objectives, Gori-Giorgi et al.|3a| 
have already proposed this formula for determining an 
average density for an inhomogeneous density, in the 
context of pair-distributions in 2-electron atoms. We 
have also shown the utility of the CDA in estimating 
the correlation energy of the 2DES in the Wigner-crystal 
phase [33]. The w of the CDA is somewhat different from 
the "thickness" 3/6 often assigned to the FH distribu- 
tion. In fact, the constant-density slab width w for the 
Fang-Howard b is given by w = 16/(36). The quasi-2D 
potential for a CDM of width w is given by 

W(r) = V(r)F(s), s = r/w, t = ^/(l + s 2 ) (13) 
1-t 



F(s) = 2s 



log- 



1 - t 



(14) 



W(r) = 



dz\dzin(z\)n{z2) 
[r 2 + (zi-z 2 ) 2 ]i/2 



(10) 



This potential tends to 1/r for large r, and behaves as 

2 / 2w r 
- In — + 1 

w \ r w 

for r < w. Thus the short-range behaviour is logarithmic 
and weaker than the Coulomb potential. The fc-space 
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FIG. 2: (a)The Fang-Howard density and the equivalent 
constant-density slab, (b) Comparison of the 1/r form with 
the quasi-2D potential from the Fang-Howard distribution for 
the FH parameter 6 = 0.1, and 0.3 and the equivalent slabs 
at w — 53.3 and 17.8 a.u. respectively. 



form of the CDM potential is: 

V usm (k,w) = V(k)F(p), p = kw (15) 
F(p) = (2/p){(e-P-l)/p+l} (16) 

The form factors F (s) and F(p) tend to unity as w — > 0. 
These r-space and fc-space analytic forms of the CDM 
lead to analytic formulae for the FH form. In our work 
we assume that a given distribution has been replaced 
by an equivalent uniform-slab distribution, and only the 
final W(r) potential enters into the exchange-correlation 
calculations (the numerical work has been checked via di- 
rect calculations as well). In the case of GaAs-HIGFETS, 
if rid could be neglected, the r s parameter specifies the 
b parameter and hence the width w of the CDM. Then 
b ~ r7 2/3 and w = 2.09494r 2/3 . 



EXCHANGE FREE ENERGY FOR QUASI-2D 
LAYERS. 



The exact exchange free energy F x for 2D layers of fi- 
nite thickness can be readily evaluated using the quasi-2D 
potential W{r) and the noninteracting pair-distribution 
functions g®j{r) of the 2-D fluid. Only the parallel-spin 



case i = j is relevant. Also, g®j{r) for a slab of finite 
thickness is identical to that for an ideally thin 2D layer, 
both at T = and at finitc-T. In fact, we find that the T 
dependence of the F x of layers of finite thickness is very 
close to that of the ideally thin case. 



Ideally-thin layer 

The first-order unscreened exchange free energy F x 
consists of F*, where i denotes the two spin species. At 
T = these reduce to the exchange energies: 



Ef/n: 



1/2 



3<A 



(17) 



Here m = n(l + C)/2, and n 2 = n(l - Q/2, Then the 
exchange energy per particle at T = 0, i.e., the total 
E x /n becomes 



E a! /n=(E?+E%)/n 



3nr s 



3/2 , 3/2, 



(18) 



where ci and ci are the fractional compositions (l±£)/2 
of the two spin species. 

We define a reduced temperature t — T/Ep, Ep — nn, 
and the species-dependent reduced chemical potentials 
fJ%/T by rji, reduced temperatures t\ = t/(l + () and 
t 2 = i/(l — C)i based on the two Fermi energies Ep\ and 
Ep2 which are Ep(l ± Q. Then we have: 



Ff/Ef = ^tf 2 



m I 2 1/2 (u)du 



, (m - «) 1/2 

The I-i/2 is the Fermi integral defined as usual: 

dxx v 



Iu(z) = 



The rji are given by 



rji = log(e 



1 



i/ti 



1) 



(19) 



(20) 



(21) 



In the paramagnetic case Eq. ED reduces to the result 
given by Isihara et al. (see their Eq.3.4; they use a dif- 
ferent definition of the Fermi integral). For small values 
of t, the exchange energy is of the form, 

E x (r s ,t) = E x (r s ,0)[l + (n 2 /16)t 2 log(i) - 0.56736< 2 + . . .] 

(22) 

The total exchange free energy is F x = Si 7 ^. The accu- 
rate numerical evaluation of Eq. 1191 requires the removal 
of the square-root singularity by adding and subtracting, 
e.g, I 2 (—\r]\)/(v — H) 1 / 2 for the case where r\ is negative, 
and v = u, and so on. 

A real-space formulation of F x = Ff + F% using the 
zeroth-order PDFs fits naturally with the approach of 
our study. Thus 



F x /n 



2irrdr 



i<j 



(23) 
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Here h%{r) = <?? 



1. In the non-interacting system 



at temperature T, the antiparallel h® 2 , viz., 5i 2 (^, T) — 1, 
is zero while 



1 



^?i( r ) = 2 Sk i. k 2 rl (' c i) n ( fc 2)e 



:(ki-k 2 )-r _ 



= -[/Ml' 



Here k, r are 2-D vectors and n(fc) is the Fermi oc- 
cupation number at the temperature T. At T = 
/(r) = 2J 1 (fc i r)/£; I . where Ji(x) is a Bessel function. As 
a numerical check, we have evaluated the exchange free 
energy by both methods, i.e., via k-space and r-space for- 
mulations. 

We present a convenient analytic fit to the exchange 
free energy which is a universal function F x (t)/E x , for 
arbitrary £. That is, the same function applies to any 
component, on using the reduced Fermi temperature of 
the spin species. The total F x is obtained by adding both 
spin contributions. The analytic fit is: 



ff(t, 0/^(0 = 

1 + Citi + C 2 t? 

1 + c 3 ti + c 4 t? 



(24) 



tanh(l/ v / <i 



The fit coefficients C t are 3.27603, 4.81484, 3.33100, 
6.51436. The temperature t{ is t/(l ± (), appropriate 
to the spin polarization. The exchange effects in the 
2DES decay more slowly with temperature than in the 
3D case where a tanh(l/t) factor appears in Eq. 3.2 of 
ref. |3 The above form does not explicitly contain the 
low-temperature logarithmic term, but it reproduces the 
value of 0.99382 at t = 0.05, while the numerical integra- 
tion gives 0.9939497. Similarly, at t = 1, 10 and 30 the fit 
(integral) returns 0.63839 (0.63839), 0.22999 (0.22990), 
and 0.13421 (0.13410) respectively. 



Thick 2-D layers. 

The T = exchange energy is modified by the layer 
thickness w. The expression for g®j(r) depends only on 
x = r/r s . Similarly, the quasi-potential W(r) depends on 
the reduced variable s — r/w. Hence the exchange energy 
of the quasi-2D layer depends only on the ratio w s — 
w/r s . The exchange energy per electron at a density n, 
given by r s , polarization £, in a layer of width w is given 
by: 

i r°° 

Ei x (r s ,Q,w) = -nr s 2ndxh°(x, QF(x/w s ) (25) 



Even though we have analytic forms for W(r) and h°(r) 
as well as their Fourier transforms, we have not found 
a convenient analytic result for the exchange energy at 
T = 0. However, the results can be parametrized by 
simple analytic-fit formulae: 

E x (r s ,C,w) = E x (r s ,(,0)Q(w s ,(), w s = w/r s 



TABLE I: Parameters fitting the exchange energy ratio as a 
function of the layer width ratio w s = w/r s . 



c 


A 


B 


G 


0.0 


0.155294 


0.142486 


0.320735 


0.5 


0.184536 


0.169822 


0.370001 


1.0 


0.19838 


0.179455 


0.462789 



Q(wsX) 



l + A(C)Vw s 



l+B((Ww s + C{()w s 



Here E x (r s , £, 0) is the exchange energy of the ideally thin 
system given by ea. 1181 The ratio Q{w s , £) is a measure of 
the reduction in the exchange energy due to the thickness 
effect. Since the effect depends on w s — w/r s , for a given 
thickness, the effect is greater for high density samples. If 
the depletion density rid in HIGFETS could be neglected, 
and if the exchange-correlation energy E xc is not included 
in the energy minimization which determines the Fang- 



Howard parameter b, then w s ~ 2.09r. 



-1/3 



The inclu- 



sion of E xc in self-consistcntly determining b changes b by 
~ 2% for low r s , but the effect becomes less important at 
higher r s . At r s =l, and 10 for £ = 0, the ratio Q is 0.652 
and 0.794 respectively. The reduction from the ideally- 
thin 2D form is clearly substantial. The exchange free 
energy F x (r s , £,T,w) at finite-T, for layers with thick- 
ness w is found to be adequately approximated by the 
temperature factor of the ideally-thin system. However, 
in calculating the effective mass from the finite-T ener- 
gies, we make independent calculations near T = and 
do not use the fit given here. 

The exchange energy for a HIGFET with nd = is 
shown in Fig. as a function of r s for £ = and temper- 
ature T/E F =0 and 0.2. 



THE EXCHANGE-CORRELATION ENERGY 
FOR QUASI-2D LAYERS. 

The correlation function h°(r) yields exact exchange 
energies for arbitrary layer thicknesses. In contrast, the 
correlation energies require a coupling-constant integra- 
tion over the pair functions g{r, £, w, A) calculated using 
the quasi-2D potential XW(r) for each A. These g(r, £, w) 
can be calculated using the CHNC. On the other hand, 
the unperturbed- g approximation, found to be useful in 
Quantum Hall effect studies 0, has been exploited by 
De Palo et al.[l2|. They have used the pair functions 
g(r,£,w — 0) of the ideally thin layer obtained form 
QMC to calculate a correction energy A given by, 

AE = (n/2) J 2irrdr[W(r) - V(r)]h(r, (,w = 0) (26) 

Then the total exchange-correlation energy E xc (r s , £, w) 
is obtained by adding to A the known E xc of the ideally- 
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FIG. 3: (a)The exchange energy E x (Hartree a.u.) of a 
2DES in a HIGFET compared to that of an ideal 2DES, at 
T/E F = and 0.2, with ( = 0. (b)The correlation energy E c 
at T = and ( = 0, for the HIGFET layer. HIGFET(upg) is 
the "unperturbed-5" approximation. HIGFET(CHNC) is the 
full calculation. This is compared with the correlation energy 
of a 3d slab model, and the "slab+rod" model. The QMC 
datum for a HIGFET ar r s = 5 is from De Palo et al[3^|. 



thin system. The above equation can also be applied 
at finite temperatures using the finite-T pair functions 
g(r s ,£,T) obtainable from the CHNC procedure. 

De Palo et al^3| have performed Diffusion-Monte- 
Carlo simulations at r s =5 for HIGFETS with b = 
0.8707, i.e, a CDM width iu=6.1256 a.u., and find that 
the error in this approach compared to the full simula- 
tion is about 2%. This full QMC result at r s = 5 is 
shown in the lower panel of Fig. [3J Since the HIGFET 
system approximates to a thin-layer as r s increases, this 
approach is probably satisfactory for r s > 5. The method 
becomes unreliable for small r s , and definitely fails below 
r s = 2. Also, the "unperturbed-g" approximation fails to 
include the renormalization of the kinetic energy picked 
up via the coupling constant integration over the fully 
consistent g(r,w). We report results (Fig. I2J) from the 
full coupling-constant integration of the g(r,w), (Fig. 
lower panel, CHNC) as well as from the "unperturbed-g" 
approximation [1 2} used by De Palo et al. 

In parametrizing the quasi-2D correlation energy 
E c {r s ,w), we present an intuitive model of E c (r s ,w). 
For small r s , the ratio w/r s is large and the electrons 



are like 1-D wires with the axis normal to the 2D plane 
and interacting with a log(io/r) interaction (c.f., Ea. ll3fl . 
However, at large r s we have 3-D like electron disks with 
w and r s of comparable magnitude in the density regimes 
of interest in HIGFETS. Thus we model the quasi-2D E c 
as an interpolation between a ID like form and a 3D like 
form. First we consider a purely 3D model. Given the 
2D-density r s and an effective CDM width w, we define 
an effective 3D density parameter , purely for calcu- 
lating its correlation energy. It will be seen that this 3D 
model is excellent for r s > 7. When r s becomes small (i.e, 
less than ~ 3), the effective width of the 2D layer, viz., 
w/r s becomes large and a ID log-interaction model j^l 
is needed. To capture the rod-like regime, we define the 
"rod like" correction AE C for r s < 7 by: 



AE C (( = 0) = a a + a L log(l/r s ) + a x r s 



(27) 



where, for HIGFETS, a L = 0.0221788, a x = 0.00365169, 
and a = 0.0192979 for C = 0. This AE C is added to the 
3-D slab form given below. The fit parameters for the ( = 
1 "rod-like" correction are a =0.013337, a L =0.0084787, 
and ai=0.0006821, to be applied for r s < 15. 

For the 3D slab-like regime (i.e, r s > 7 for C=0, r s > 15 
for £ = 1) we define a (^-dependent 3-D density parameter 
and a correlation energy via: 



E c (r s ,C,w) 

r 3d 
' s 

z(0 



E 3 c d (r 3 s d ,C) 
r s /F{r s ) 

2^2 



(4-C)Va 



1.5 



+ d- 5 



(28) 
(29) 
(30) 

(31) 



The 3D correlation energy Ef d (rl d , Q is that given 
by, e.g., Ceperley and Alder |38j, or Gori-Giorgi and 
Perdew 39] . We see from the lower panel of Fig. |3 that 
the correlation energy for small r s , calculated using the 
3D slab begins to go below the "unperturbed-g" approx- 
imation of Eq. 1261 consistent with the trend of the ideal 
2D gas, and the trend of the only QMC data point avail- 
able for a HIGFET, at r s — 5. The exchange-correlation 
energy obtained from the full CHNC calculation is in 
excellent agreement with the QMC datum. The curve 
labeled "Slab+rods" in Fig. O is the combined formula, 
Eq. 123 and Eq. |2H1 for the correlation energy/electron, 
with, for example, the region r s ~ 7 for £ = obtained 
by linear interpolation between r s =6 and r s = 8. This 
is clearly seen to reproduce the full CHNC results very 
well. 



Correlation energy at finite temperatures. 

The correlation contribution to the Hclmholtz free en- 
ergy of the ideal 2-D layer, and layers of thickness w can 
be easily calculated using the approach of Eq. where 
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TABLE II: Low-temperature data for the exchange- 
correlation contribution to the Helmholtz free energy, F xc per 
electron in atomic units at r s — 5, for the ideal 2D system and 
for a CDM of width 6.1256 au., i.e., Fang-Howard 6=0.8707. 
The packing fractioin r\ which defines the hard-disk bridge 
function is also given 



t = T/Ef 


V 


F xc (ideal) 




F IC (CHNC) 


0.05 


0.3718 


-0.16902 


-0.11197 


-0.11467 


0.08 


0.3725 


-0.16883 


-0.11177 


-0.11448 


0.10 


0.3730 


-0.16867 


-0.11161 


-0.11433 


0.12 


0.3733 


-0.16850 


-0.11144 


-0.11417 



the CHNC-generated finite- T pair functions are used. A 
typical set of results at very low temperatures is given in 
Table [H] Here we have also given the packing fraction 
r] of the hard-sphere bridge function used to mimic the 
three-body and multi-body correlation contributions. As 
discussed in earlier work 27] , the F x and the F c at very 
low T contain logarithmic terms which cancel with each 
other, so that the sum F xc — F x +F c is free of such terms. 
From our numerical work we find that the T dependence 
of the F x and F c of layers of finite thickness is very sim- 
ilar to that of ideally thin layers. Hence we assume that 
the logarithmic corrections are also similar. At r s = 5 
the cancellation is good to about 75%, and this improves 
as r s increases. Although the two-component fluid (up 
spins and down spins) involves three distribution func- 
tions, we have, as before j^], used only one hard-disk 
bridge function, -B12, as clustering effects in ga are sup- 
pressed by the Pauli-exclusion. However, at high densi- 
ties (low r s ), the use of three bridge functions seems to 
be needed for satisfying various subtle features that are 
needed to ensure the exact cancellation of logarithmic 
energy terms etc. Instead of introducing additional fea- 
tures into the CHNC method, we have however retained 
the single bridge-function model that was used by us so 
far [12. 



r 

s 




r 

s 



FIG. 4: The Energy difference between the fully polarized and 
unpolarized states of 2D layers with different CDM widths w. 
The HIGFET w varies as r 2 J^ and is ~ 18 near the spin-phase 
transition. 

w is only 18.4, and the SPT remains intact and occurs 
at a somewhat higher r s , as shown by de Palo et alflij]. 
and also in Fig. 0] A natural consequence of delaying the 
SPT is to decrease the spin-susceptibility enhancement. 
The effective thickness of the quasi-2D layer can be in- 
creased by suitably designing the shape of the potential 
well, or including an additional subband, and in this case 
the SPT can be circumvented. However, a discussion of 
higher subband effects is outside the scope of this study. 



The transition to a spin-polarized phase. 

QMC simulations as well as CHNC calculations show 
that there is a spin-polarization transition (SPT) in the 
ideally-thin 2D electron fluid near r s ~ 26. On the other 
hand, the correlation contributions dominate over the ex- 
change energy in the 2-valley2D system in Si MOSFETS, 
and the SPT is suppressed [25j. The rapid increase in to*, 
with g* remaining unchanged while r s is increased, ob- 
served by Shashkin et al.|40j was found to be consistent 
with this picture jl^]. In finite-thickness 2D layers, as 
the CDM width w increases, the location of the SPT is 
pushed to higher values, as seen in Fig. QJ In the case 
of HIGFETS used by, e.g., Tan et al.|nJ7the width w 
increases with r s , but only as r^ 3 . Thus at r s =26, the 



THE SPIN-SUSCEPTIBILITY, EFFECTIVE 
MASS AND THE ^-FACTOR 

The results for the exchange-correlation free energy 
F xc {rs : C, T) for the ideal 2D system and the thick-layer 
system contain all the information needed to calculate 
the spin-susceptibility enhancement, the effective mass 
to* and the effective Lande factor g* , for the ideal sys- 
tem and the thick layer. In fact, the following quantities 
are calculated from the respective second derivatives of 
the energy. 



[d 2 F xc (t)/dt 2 
[d 2 F Q (t)/dt^ 



(32) 
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□ HIGFET, this work. 
— HIGFET, Zhu et al. 
Ideal-2D, QMC 




Ideal 2D Layer 



x — x PT, Zhang et al. 

CHNC, this work 

/ g- -o G + &G /D, Asgari et al. 

□ QMC, Kwon et al. 




FIG. 5: (a)The spin-susceptibility enhancement, i.e., 
Xs/X-P— m *9* as a function of r s for the HIGFET. The ex- 
perimental data of Zhu et al.Q are also shown. (b)Variation 
of m*g* as a function of the spin-polarization for a HIGFET. 



Xp/Xs 



(m g 



[d 2 F xc {Q/de 

[d^F (C)/dC 2 ] 



(33) 



We use the available QMC results for the ideal 2D 
exchange-correlation energy at T = 0, and where con- 
venient, the QMC pair-distribution functions at T = 
as parametrized by Giri-Giorgi et al[3(|. The CHNC is 
used to obtain the pair-functions for situations (e.g, at 
finite-T and finite thickness) where the QMC data are 
simply not available or difficult to use. In most cases, 
replacing the QMC-PDFs with the CHNC ones, or using 
the "unperturbed-g" approximation leads to relatively 
small changes. The exception is in the calculation of to* , 
where the "unperturbed-g" approximation, Eq. 1261 is in- 
adequate. 



The effective mass m*. 

In Fig. HJ1 we present our results for the ideal-2DES. 
No experimental results are available for this case, but 
limited QMC simulations 01 as we ^ as results from di- 
agrammatic perturbation theory (PT) 0, ^| are avail- 
able. The CHNC based to* values have an error of at 
most ± 2%. The ideal 2D-layer to* shows a rapid rise 
around r s =2 to 5, in good agreement with the four val- 



FIG. 6: The effective mass m* of an ideal 2D layer (w = 0) 
obtained from CHNC are compared with the the Quantum 
Monte Carlo data of Ref. Il7ll and the perturbation theory 
calculations of Zhang et al.|l4ll. and Asgari et al. |13jl . i.e., 
their G+&G-/D calculation. 



ues from QMC, and then slows down in strong contrast 
to the to* proposed by Asgari et al., and Zhang et al. We 
have displayed the model denoted G+&G_/ D by Asgari 
et al., as being their optimal choice from among the many 
models given in Ref. 0], where they also contest the 
analysis of Zhang et al. The PT approaches are partly 
scmi-phenomenological in that QMC data are input into 
local-field factors and other ingredients of the calculation; 
the choice of the vertex functions, treatment of on-shell 
or off-shell effects, whether to use self-consistent propa- 
gators or not, etc., are components of the prescription 
used by different workers. The strong disagreement of 
the PT-based to* with the QMC-based to* is notable. 

The CHNC method has some similar uncertainties, es- 
pecially in the use of a Percus-Yevik hard-sphere bridge 
function B(r) to capture the back- flow like three-body 
contributions to the PDFs and the total energy. As seen 
in Fig. ^ the T dependence of B(r) seems quite small, 
and our initial calculations of m*, reported in Ref. |42l | 
were based on the zero-T form of B(r). This leads to an 
to* which drops slightly below unity and remains there. 
In this calculation we have used the proper T-dependent 
bridge function and the calculated m* is in good agree- 
ment with the QMC-based m* . This might be somewhat 
coincidental, as the QMC results are also based on sensi- 
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FIG. 7: The effective mass m* H in a HIGFET. The PT results 
of Zhang et al., and Asgari et al., for the HIGFET are close 
to the ideal-2D results from QMC. The "unperturbed-g" ap- 
proach to the energy using CHNC-based PDFs gives a slight 
reduction from the ideal 2D. For the full CHNC calculation, 
see Fig. 



tive approximations. However, it means that we do have 
a B(r) which is consistent with current QMC results, 
and hence may be used with greater confidence in study- 
ing thick-2D systems. Another point in favour of our 
model of B(r) is seen in the local-field factor (LFF) of 
the ideal 2DES response function. A study of the LFF of 
the 2D response |2fj shows that the formation of singlet- 
pair correlations is essentially complete by r s ~ 5, and 
after that the structure of the fluid remains more or less 
unchanged, until the SPT is reached. The hard-sphere 
model of B{r) provided a satisfactory description of the 
short-ranged features of the 2DES-LFF. The rapid rise 
in to* up to r s ~ 5 and the subsequent slow-down is 
probably related to the formation and persistence of the 
singlet structure in the 2D fluid revealed by the form of 
the LFF. 

In Fig. we present a comparison of various theoreti- 
cal calculations of the effective mass m* H of the electrons 
in the HIGFET. The PT calculations of Zhang et al, 
and Asgari et al., show a strong decrease of m* from 
the PT values in the ideal 2DES. Our calculations, us- 
ing the "unperturbed-g" approximation, Eq.ESl lead to a 
m* H which is only slightly reduced by the thickness effect. 
This m* H curve is in close to that of Asgari et al. This is 



clearly a numerical accident. According to Asgari et. al., 
the difference between the ideal to* and m* H increase as r s 
increases. In our calculation using the "unperturbed-i?" 
approximation, the difference, already quite small, seems 
to diminish as r s increases. In fact, as r s increases, the ra- 
tio w/r s of the HIGFET layer decreases and the thickness 
effect may be expected to decrease, unless the difference 
between to* and m* H is driven by some other effect. This 
"other effect" is revealed by giving up the "unperturbed- 
g" approximation, and using the full thick-layer 2DES 
pair-distribution function at finite T, calculated using the 
CHNC, to evaluate the total free energy F(r B ,T) of the 
quasi-2D system, and hence the m* H . In Fig.|Hlwe display 
the PDF of the quasi-2DES of a HIGFET at r s = 5, and 
compare it with the PDF of the ideal 2DES. The differ- 
ence between the ideal and quasi systems is embodied in 
the form factor F(r). The reduced Coulomb repulsion 
at small-r leads to a large pile-up of electrons around the 
electron at the origin. This means, the electron has to 
drag this charge pileup and this contributes an enhanced 
to* to the thermodynamic and transport properties of 
the quasi 2DES. 

The experimental results of Tan et. al., for to* show 
an increase of ~ 150% between r s =3 and r s = 6. Our 
results from the full CHNC calculation, as well as the 
perturbation results of Asgari et al., are shown in the 
lower panel of Fig. |H1 Given the failure of the PT calcu- 
lations to reproduce the QMC-based to* for the ideal 2D, 
it is difficult to evaluate the reliability of the PT-based 
m* H . The PT-overestimate of to* of the ideal 2D is typi- 
cal of the RPA-like character of these theories which are 
likely to predict spin transitions at relatively high densi- 
ties. Also, we believe that if the same PT prescriptions 
were used to evaluate the one-top value g(0) of the PDFS 
of the 2DES and the quasi-2DES, another measure of the 
short-comings of the PT methods would be revealed. 

Enhanced spin susceptibility and the Lande-g factor. 

The product m*g* is given by the ratio of the static 
spin susceptibility \s to the ideal (Pauli) spin susceptibil- 
ity xp- The long wavelength limit of the static response 
functions are connected with the compressibility or the 
spin-stiffness via the second derivative of the total energy 
with respect to the density or the spin polarization|25j. 
De Palo et al.0 have calculated m*g* from the QMC 
pair distribution functions and shown that they obtain 
quantitative agreement with the data for very narrow 2- 
D systems |4l|as well as for the thicker systems found 
in HIGFETSQ. The CHNC PDFs are close approxi- 
mations to QMC results, and when used in Eq. [23 yield 
correction energies which are in goo d ag reement with the 
energies obtained by De Palo et al|35l|. In Ref. we 
showed that the rapid enhancement of m*g* in Si/Si02 
2DES is a consequence of the increase in to* with r s , 
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FIG. 8: (a)The spin-averaged g(r) at r a = 5 and T/E F =Q.l, 
for an ideal 2DES, and for a HIGFET calculated using CHNC 
(b)The effective mass m* H in a HIGFET from fully self- 
consistent CHNC calculations finite-T. The experimental 
data are from Tan et al^l. 

and that the g* does not increase with r s because there 
is no spin-phase transition in the 2-valley case. In the 
HIGFET system there is a slightly delayed SPT, as seen 
in Fig. [S] Hence m*g* increases with r s , while to* also 
increases quite rapidly, due to the enhanced "on-top" cor- 
relations shown (Fig. [SJ in the PDF of the quasi-2DES. 
The resulting g * of the HIGFET is shown in Fig. 



CONCLUSION. 

We have presented a detailed study of the effect of 
many-body interactions in quasi-2D electron layers using 
a single theoretical framework which involves the calcu- 
lation of the pair-distribution functions of the system via 
a classical representation of the quantum fluid. A proce- 
dure for replacing the inhomogeneous transverse distri- 
butions via a constant-density model, i.e., an equivalent 
homogeneous distribution, has also been demonstrated. 
Easy to use parametrized fit formulae for the exchange 
energy at zero and finite-T have been presented. A sim- 
ple numerical scheme for calculating the correlation en- 




FIG. 9: The ideal 2DES g-factor is obtained from the QMC 
m*g* of ref. Q, divided by the CHNC ideal 2D m*. The 
experimental HIGFET t/Jj (boxes) is from the m*g* of Zhu 
et al., divided by the m* of Tan et al. The inset shows the 
spin-susceptibility enhancement m*g* from the E xc (r s ,Q cal- 
culated from Eg. 1261 where the ideal 2D g(r,w — 0) is used, 
and from the full CHNC calculation using the g(r,w) consis- 
tent with the quasi-2D potential. 



ergy of a thick 2D layer,, via a 3-D slab model combined 
with a 1-D rod model, has also been demonstrated. We 
find that the thickness effect on the spin-phase transi- 
tion etc., provides a clear picture of the changes in the 
spin-susceptibility enhancement leading to a strong in- 
crease in the g-factor, while to* is increased due to the 
enhancement of the "on-top" correlations arising from 
the reduction of the Coulomb potential in thick layers. 
However, unlike in the case of the effective mass data 
for Si/Si02 systems [27tl4fj| . these results do not provide 
good quantitative agreement with the effective-mass data 
for HIGFETS recently reported by Tan et al. This may 
be due to our use of the ideal 2DES bridge function even 
for the HIGFETS. 
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